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ABSTRACT 

We use semi-analytic models implemented in the Millennium Simulation to analyze 
the merging histories of dark matter haloes and of the galaxies that reside in them. 
We assume that supermassive black holes only exist in galaxies that have experienced 
at least one major merger. Only a few percent of galaxies with stellar masses less than 
M* < 10 10 M Q are predicted to have experienced a major merger and to contain a black 
hole. The fraction of galaxies with black holes increases very steeply at larger stellar 
masses. This agrees well with the observed strong mass dependence of the fraction of 
nearby galaxies that contain either low-luminosity (LINER-type) or higher-luminosity 
(Seyfert or composite-type) AGN. We then investigate when the major mergers that 
first create the black holes are predicted to occur. High mass galaxies are predicted 
to have formed their black holes at very early epochs. The majority of low mass 
galaxies never experience a major merger and hence do not contain a black hole, but a 
significant fraction of the supermassive black holes that do exist in low mass galaxies 
are predicted to have formed recently. 
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1 INTRODUCTION 

By studying active galactic nuclei (AGN), we learn about 
the physical mechanisms that trigger accretion onto the cen- 
tral supermassive black holes of galaxies. When a black hole 
accretes, it increases in mass. By studying populations of 
AGN at low and at high redshifts, we hope to infer the his- 
1 tory of how black holes build up their mass. 

It has been established that s upermassive black holes 
most occur in galaxies with bulges l|Kormendv fc Richstonel 

1 19951 ). and that the mass of the black hole correlates with 
the lumin osity and the stellar velocity dispersion of the 
host bulge dMagorrian et al .1 19981 : iFerrarese fc Merrittj200d ; 

iGebhardt et al.l l2000h . This indicates that the formation 
of galaxies and supermassive black holes are likely to be 
closely linked. In the local Universe, the fraction of bulge- 
dominated galaxies hostin g AGN decreases at low er stel- 
lar masses i|Ho et al.l fl997l : iKauffmann et alJlioOol ). In or- 
der to form a black hole, it is necessary for gas to lose 
angular momentum and sink to the centre of the galaxy 
l|Haehnelt fc Reeslll993l ; IVolonteri et al.ll2003l ). The gravita- 
tional torques that operate during galaxy-galaxy mergers 
are known to be a very effective mechanism for concentrat - 
ing gas at the centers of galaxies (|Mihos fc H crnquist 19961 ). 
Models for AGN evolution have often assumed that black 
holes are formed and fuelled, and AGN activity is triggered 



Email: wanglan@mpa-garching.mpg.de 



during major mergers of galaxies (IKauffmann fc Haehneltl 



during major mergers ot galaxies (Kaurmiaim 
l200d : IWvithe fc Loeblliool ICroton et alB ooe) 



At low and moderate redshifts, there is no conclusive 
observational evidence that mergers play a significant role 
in trigge r ing A GN activity in galaxies. In the local Universe, 
iLi et al.l l|2006l ) have shown that narrow line AGN do not 
have more close companions than matched samples of inac- 
tive galaxies. Even at intermediate redshifts (z ~ 0.4 — 1.3), 
moderate luminosity AGN hosts do not have morphologie s 
indicative of an ongoing merger or interaction l|Hasanll2007l ). 
The conclusion seems to be that although major mergers 
may be responsible for AGN activity in some galaxies, other 
fueling mechanisms are likely to be most important in the 
low redshift Universe. It has also been established that high 
mass black holes have largely stopped growing at early cos- 
mic epochs, whereas low mas s black holes are still accret- 
ing at significant rates today (|Heckman et al.ll2004l ). X-ray 
observations show that very high-luminosity AGN activity 
peaked at early cosmic epochs {z ~ 2), while low-luminosity 
AGN activity peaks at lower redshifts (jSteffen et al.l 120031 ; 
iBarger et aij|2005l : lHasinger et~aT1l2005l ). 

It has been postulated that this so-called "anti- 
hierarchical" growth of supermassive black holes can be 
explained if there are two modes of accret ion onto black 
holes that have very di fferent efficiencies l|Merlonil 120041 : 
I Mueller "fc~ Hasingcr 2007f) . The early formation formation of 
"new" black holes may result in very luminous quasar-like 
events. To form a supermassive black hole, a more violent 
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process such as a major merger may be required to funnel 
a large amount of gas into the central region of the galaxy. 
Subsequent accretion of gas onto already existing black holes 
may b e an inefficient process and produce lower luminosit y 
AGN (|Haehnelt fc Reeslll993l ; lDuschl fc Strittmatterll2002l ). 
The history of accretion after the black hole is formed may 
not necessarily be tightly linked to the dynamical history 
of the galaxy, but may be controlled by the accretion and 
feedback processes occurring in the vicinity of the black hole 
itself. 

In this work, we use the combination of the Millennium 
Simulation and semi-analytic models of galaxy formation to 
study the fraction of galaxies that have undergone major 
mergers as a function of mass and cosmic epoch. We inves- 
tigate whether this can be related to the demographics of 
black holes in the local Universe and to the apparent dis- 
appearance of the most luminous quasar activity in massive 
galaxies at late times. 

In Sec. [51 we briefly introduce the simulation we use 
and explain how galaxy mergers are tracked in the simula- 
tion. In Sec. [3] we show that if we assume that black holes 
only form when galaxies undergo major merging events, then 
most present-day low mass galaxies are predicted not not to 
contain black holes and hence will not host AGN. In Sec. |4j 
we use the simulations to predict when galaxies of different 
masses have underdone their first major merger. Conclusions 
and discussions are presented in the final section. 



2 SIMULATION AND MERGER TREES 

The Millennium Simulation (|Springel et ai1 l2005) is used in 
this work to study the merging histories of dark matter 
haloes. The merging histories of galaxies can be inferred 
when the simulation is combined with semi-analytic models 
that follow gas cooling, star formation, supernova and AGN 
feedback and other physical processes that regulate how the 
baryons condense into galaxies. 

The Millennium Simulation follows N = 2160 3 particles 
of mass 8.6 x 10 8 /i _1 Mg from redshift z — 127 to the present 
day, within a comoving box of 500/i _1 Mpc on a side. The 
cosmological parameters values in the simulation are consis- 
tent with th e determinations fr om a combined analysis of 
the 2dFGRS(IColless et alj|200 ll ) and first year WMAP data 
(|Spergel et al.ll2003D . A flat ACDM cosmology is assumed 
with fi m = 0.25, fib = 0.045, h = 0.73, fi A = 0.75, n = 1, 
and as = 0.9. 

Full particle data are stored at 64 output times. For 
each output, haloes are identified using a friends-of-friends 
(FOF) group-finder. Substructures (or subhaloes) within a 
FOF halo are locat ed using the SUBFIND algorithm of 
ISpringel et"aH ^OOll ). The self-bound part of the FOF group 
itself also appear in the substructure list. This main subhalo 
typically contains 90 percent of the mass of the FOF group. 
After finding all substructures in all the output snapshots, 
subhalo merging trees are built that describe in detail how 
these systems merge and grow as the universe evolves. Since 
structures merge hierarchically in CDM universes, for any 
given subhalo, there can be several progenitors, but in gen- 
eral each subhalo only has one descendant. Merger trees are 
thus constructed by defining a unique descendant for each 
subhalo. We refer below halo to the main substructure that 



can represent the FOF halo, while subhalo refers to substruc- 
ture other than the main one. Halo merger happens when 
two FOF group merge into one group and one of the haloes 
becomes a subhalo of the larger structure. 

The substructure merger trees form the basic in- 
put to the semi-analytic model used to associate galax- 
ies with haloes/subhaloes l|De Lucia fc Blaizotl [2007). The 
semi-analytic galaxy catalogue we are using in this study 
is publicly available. A description of the publicly available 
catalogues, and a link to the database can be found at the 
webpage: http://www.mpa-garching.mpg.de/miilennium7] 
Once a halo appears in the simulation, a (central) galaxy 
begins to form within it. The central galaxy is located at the 
position of the most bound particle of the halo. As the sim- 
ulation evolves, the halo may merge with a larger structure 
and become a subhalo. The central galaxy then becomes a 
satellite galaxy in the larger structure. The galaxy's position 
and velocity are specified by the position and velocity of the 
most bound particle of its host halo/subhalo. Even if the 
subhalo hosting the galaxy is tidally disrupted, the position 
and velocity of the galaxy is still traced through this most 
bound particle. Galaxies thus only disappear from the simu- 
lation if they merge with another galaxy. The time taken for 
a galaxy without subhalo to merge with the central object 
is given by the time taken for dynamical friction to erode its 
orbit, causing it to spiral into the centre and merge. This is 
calculated using the standard Chandrasekhar formula. All 
the information about the formation and merging history of 
galaxies is stored. 

By analyzing these halo and galaxy merger trees, we are 
able to track when two haloes merge together and whether 
the galaxies within them also merge into a single object by 
the present day. In this study, we focus on mergers between 
satellite and central galaxies, and excl ude mergers between 
two satellites. These events are rare (ISpringel et al.ll200ll ) 
and neglecting them should not affect our conclusions about 
the incidence and fueling of black holes in galaxies. 



3 HALO AND GALAXY MERGERS 

In this study, we assume that black holes form when a galaxy 
undergoes a major merging event. Galaxies that have never 
experienced a major merger do not have a black hole. We 
define major mergers as events in which the mass ratio of 
the two progenitors is greater than 0.3. For halo merger, the 
mass ratio is the virial mass ratio of two progenitor haloes. 
For galaxy merger, it is the stellar mass ratio of two progen- 
itor galaxies. When we track mergers in the simulation, we 
include major mergers that occur in all branches of the tree, 
not just the "main branch" . 

Since galaxies reside in dark matter haloes and are able 
to merge only once their host haloes have coalesced, we begin 
by analyzing the merging histories of the dark matter haloes 
themselves. In the left panel of Fig. [TJ we plot the average 
number of major mergers a present day dark matter halo has 
experienced over its lifetime as a function of halo mass. Note 
that in this analysis we track mergers down to an effective 
resolution limit of 20 particles, which corresponds to a halo 
of mass 1.7 x 10 10 /i _1 Mq. We see that the number of major 
mergers (above the resolution limit) experienced by a halo is 
a strongly increasing function of mass; haloes with present- 
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Figure 1. Left panel: the average number of major mergers that a dark matter halo of given mass has experienced over its lifetime. 
Right panel: the fraction of haloes of given mass that have had at least one major merger. 
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Figure 2. Left panel: The relation between the stellar mass of the central galaxy and the the mass of its host dark matter halo as 
predicted by the semi-analytic models of De Lucia & Blaizot (2007). The error bars indicate the 95 percentile range in stellar mass at a 
given value of M(, a ; . Middle panel: The solid line shows the fraction of dark matter haloes that have experienced at least one major 
merger as a function of the stellar mass of the central galaxy. The dotted line shows the fraction of central galaxies of given mass that 
have had at least one major merger. Right panel: The solid line shows the average number of major mergers experienced by a dark matter 
halo as a function of the stellar mass of its central galaxy. The dashed line shows the average number of major mergers experienced by 
the central galaxy itself. 



day masses of 10 12 Mq have typically experienced only one 
one major merger, whereas the progenitors of present-day 
haloes with masses of 10 15 M Q have merged with each other 
close to 100 times. 

In the right panel, we show the fraction of haloes that 
have had at least one major merger during their lifetime, 
as a function of halo mass. The fraction of haloes that have 
had major mergers also increases rapidly with halo mass. 
Almost all haloes more massive than 10 13 /i _1 Mq have had 
at least one major merger and this fraction drops to around 
20 % for haloes with masses of around 10 11 /i" 1 M Q . 

We now investigate the fraction of galaxies that have 
had major mergers. The results are shown as a dotted line 
in the middle panel of Fig. [5] Rather than rising steeply 
as a function of mass, the galaxy major merger fraction re- 
mains close to zero up to a stellar mass of 10 10,5 Afo and 
then rises sharply. This is somewhat surprising in view of 
the behaviour of the same quantity for dark matter haloes, 
plotted in the right-hand panel of Fig. [1] For reference, we 



have plotted the relation between the stellar mass of a cen- 
tral galaxy and the mass of its host halo in the left panel 
of Fig. [2] as p redicted by the semi-analy tic models we use 
in this study (IDe Lucia fc BlaizotJ 120071 ). This mean rela- 
tion can be used to transform between central galaxy mass 
and halo mass in an approximate way (this conversion ne- 
glects scatter between the two quantities and the fact that 
some galaxies are actually satellite systems). If the fraction 
of galaxies with major mergers followed the relation derived 
for their host haloes, this would yield the solid curve in the 
middle panel of Fig. Why are the merging histories of 
galaxies and their host haloes so different? 

Once two dark matter haloes merge, the galaxies in- 
side them will merge together over a timescale that is de- 
termined by dynamical friction. Upon investigation, we find 
that nearly all galaxies that have experienced major merg- 
ers are located in dark matter haloes that have also expe- 
rienced a major merger. There are almost no galaxy major 
mergers that have occurred in a halo that has only experi- 
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Figure 3. The thick solid line (HMM) shows the fraction of cen- 
tral galaxies whose progenitor haloes have had at least one major 
merger. The other lines split this sample of central galaxies ac- 
cording to the history of the central galaxy itself. The dotted line 
(GnoM) shows the contribution from central galaxies that have 
not experienced a mergers of any kind. The dashed line (GmM) 
shows the contribution from central galaxies that have experi- 
enced only minor mergers. The thin solid line (GMM) shows the 
contribution from central galaxies that have experienced major 
majors. 



enced a minor merger (~ 0.15 percent). However, the con- 
verse is not true; we find that a substantial fraction of halo 
major mergers give rise to galaxy minor mergers. This is 
illustrated in the right-hand panel of Fig. [2] The solid line 
shows the number of major mergers experienced by the pro- 
genitor haloes of a present-day central galaxy as a function of 
its mass. The dashed line shows the number of major merg- 
ers experienced by their progenitor galaxies. As can be seen, 
the number of major mergers experienced by the progeni- 
tor galaxies is an order of magnitude smaller. Notice that 
the number of galaxy mergers is less than 1 for galaxies up 
to ~ 10 11 /i -1 Mq, and increase steeply for massive galaxies. 
This is in nice agree ment with what is shown in Fig. 9 of 
lDe~ Lucia et alj l|2006l ). which shows the number of effective 
progenitors as a function of the stellar mass for elliptical 
galaxies. 

In Fig. [3] we again plot the fraction of central galaxies 
of a given mass whose progenitor haloes have had a ma- 
jor mergers (thick solid line). The thin solid line shows the 
fraction whose progenitor galaxies have had a major merger. 
The dashed line shows the fraction of such galaxies that have 
had minor mergers and the dotted line is the fraction that 
have had no merger of any kind. The main conclusion that 
can be gleaned from this plot is that the reason why the 
thick solid and thin solid curves differ in shape, is because 
at lower stellar masses, major mergers between the progen- 



itor haloes often correspond to minor mergers between the 
progenitor galaxies. 

How can we understand this? During the period of time 
between the merger of the two haloes and the merger of the 
galaxies within them, the stellar mass of the smaller "satel- 
lite" galaxy remains about the same because ongoing star 
formation is quenched when the gas surrounding the galaxy 
is shock-heated and no longer cools onto the satellite. The 
central galaxy, however, will continue to increase in mass as 
a result of cooling and star formation. The stellar mass ratio 
of two galaxies therefore becomes smaller as a function of 
time. 

This is illustrated in Fig. [4] For every merging event 
that occurs over the history of a galaxy, we record stel- 
lar mass ratio information at the time when the progenitor 
haloes merge and at the time when the galaxies themselves 
merge together. For simplicity, we keep information for one 
randomly chosen merging event in the history of each galaxy. 
In the left panel of Fig. [4] we plot the average time that 
elapses between the time when the two haloes merged and 
the time when the galaxies themselves merged. Results are 
shown as a function of galaxy stellar mass and the the error 
bars indicate 68 percentile range in the distribution of delay 
times. As can be seen, the typical delay time is around 2 Gyr, 
but individual time delays can range between 1 and 5 Gyr. 
The delay times are typically shorter for the progenitors of 
more massive galaxies. 

In the right panel of Fig. [4] we plot the average stellar 
mass ratios of the galaxies at the time when their haloes 
merge (solid curve) and at the time when the two galaxies 
themselves merge (dashed line). Notice that the stellar mass 
ratio can sometimes be larger than 1; this happens when 
the galaxy inside the smaller halo is more massive than the 
galaxy in the larger halo. As we expect, the mass ratio of 
galaxies at the time when the galaxies merge is smaller than 
that it is at the time when the haloes merge. This effect is 
somewhat larger for the mergers that give rise to the most 
massive galaxies at the present day. 

3.1 Comparison with Observations 

In this section we have seen that the fraction of galaxies 
that have experienced one or more major merging events is 
predicted to very close to zero at stellar masses less than 
~ 1O 1O M0, but a very steeply rising fraction of stellar mass 
for M„ > 1O 1O M0. We now compare this prediction with the 
fraction of Sloan Digital Sky Survey galaxies that contain an 
AGN. We restrict the SDSS sample to redshifts z < 0.06 so 
that we are still able to de tect AGN with weak lin e emis- 
sion (LINERs). As shown bv lKauffmann et all l|200&t ). weak- 
lined AGN become progressively more difficult to identify at 
higher redshifts using SDSS spectra. This is because these 
spectra are obtained through 3 arcsecond diameter fibre 
apertures and the contribution from the stellar population 
of the host galaxy becomes increasingly dominant in more 
distant galaxies. 

The results of the comparison are shown in Fig. [S] The 
black curve shows the fraction of galaxies in the Millen- 
nium Simulation of given stellar mass that have had at least 
one major merger. The red histogram shows the fraction of 
galaxies in the SDSS survey that are classified as AGN. As 
can be seen, both fractions rise steeply from values close to 
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Figure 4. Some characteristics of the merging events with the highest stellar mass ratios that take place during the history of a galaxy: 
Left: the average time delay between the time that the progenitor haloes merge and the time that the central galaxies merge. Right: the 
stellar mass ratio of the galaxies at the time that the haloes merge (solid line) and at the time when the central galaxies merge (dashed 
line). All results are plotted as a function of the stellar mass of the central galaxy and error bars show the 68 percentile dispersion around 
the mean value. 
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Figure 5. The solid curve shows the fraction of galaxies of given 
stellar mass that are predicted to have experienced at least one 
major merger. The red histogram shows the fraction of SDSS 
galaxies with z < 0.06 that are classified as AGN. The dotted 
and dashed lines show the results from simulation when the mass 
ratio threshold for defining major merger is changed to 0.2 and 
0.4. 



zero at M» < 10 10 h~ x Mq to nearly unity at stellar masses 
greater than 10 ll /i — 1 Afg. In Fig. [S] the dotted and dashed 
lines show the results from simulation when the mass ratio 
threshold for defining major merger is changed to 0.2 and 



0.4. Compared with the solid line where we use 0.3 as the 
mass ratio threshold to define a major merger, the increas- 
ing trends are about the same for different thresholds in the 
range from 0.2 to 0.4. 



4 FIRST BLACK HOLES 

In this section, we analyze when the first major merger that 
produces the black hole in the galaxy is predicted to occur. 
In Fig. [5] we plot the distribution of the times of the first 
major merging events for galaxies with different present- 
day stellar masses. The vertical dashed lines indicate the 
median values of the distributions. The dotted bar in each 
panel indicates the fraction of galaxies in each stellar mass 
bin that have not experienced a major merger and are hence 
not included in the distribution of merging times. 

As can be seen, massive galaxies experience their first 
major merging event at earlier epochs than less massive 
galaxies. Almost no new black holes form in massive galaxies 
at the present day. The distribution of black hole formation 
times in low mass galaxies is much flatter. If the bulk of the 
black hole mass is built up in a short period following the 
first major merger, this would explain why present-day mas- 
sive black holes have stopped growing, while low mass black 
holes are still growing at a significant rate (He ckman et al.l 
l2004h . 

We now assume that the black holes formed from the 
first major mergers of galaxies can shine and be observed for 
10 7 years. By counting the numbers of such events at differ- 
ent redshifts, we can compute the evolution in the number 
density of newly formed black holes. The result is plotted 
as diamonds in Fig. [7] The comoving number density of 
such events peaks at redshift of z ~ 2 — 3, consistent with 
the observed peak in the number density of bright quasars 
l|Richards et al.l [20061 ). The decrease in the number density 
of newly formed black holes to high redshifts is less pro- 
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Figure 6. Distribution of time when a galaxy experiences its first major merger for galaxies in different stellar mass bins (solid lines). 
Galaxy stellar mass is plotted in units of h~ Mq. The vertical dashed lines show the median value of the distributions. In each panel, 
the dotted bar shows the fraction of galaxies that have never experienced a major merger. 



nounced than that found by iFan et al.l (|2001T ). who show 
that the luminous quasar density decreases by a factor of ~ 6 
from redshift 3.5 to 5. Note that we have not attempted to 
model the predicted luminosity of the quasars in this work, 
so a direct comparison with the observations is not possi- 
ble. As we have discussed, it is well possible that processes 
other than mergers contribute to the low-luminosity quasar 
population. 



5 CONCLUSIONS 

We analyze the merger histories of dark matter haloes and 
galaxies in the Millennium Simulation and use our results to 
try to understand the demographics of black holes in nearby 
galaxies. Black holes are assumed to form only if a major 
merger occurs. Although a significant fraction of low mass 
(< 10 10 Mq) galaxies have experienced minor mergers, less 
than a few percent are predicted to have experienced a major 
merger. If our assumption that a major merger is required 
in order to form a black hole is correct, the majority of 
low mass galaxies are predicted not to contain black holes 
at the present day. This is one possible exp lanation of the 
observed lack of AGN in low mass galaxies (|Ho et al.lll997l ; 
iKauffmann et alj|2003h . 

We also investigate when galaxies of different stellar 
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Figure 7. The number density of galaxies experiencing their 
first major merger (diamonds) is plotted as a function of redshift. 
Each merger is assumed to be visible for 10 7 years. 
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masses are predicted to have formed their first black holes. 
High mass galaxies form their first black holes at very early 
epochs. The distribution of formation times is almost flat 
as a function of lookback time for low mass galaxies. This 
means that if a low mass galaxy has a black hole, there 
is a significant probability that it formed in the last few 
Gigyears. We also compute the number density of newly 
formed black holes as a function of redshift. We find that the 
peak number density occurs at z ~ 2 — 3, in good agreement 
with the observed peak in the quasar space density. More de- 
tailed predictions for how AGN of different luminosities are 
expected to evolve requires a more detailed physical model 
for how the black holes accrete gas over the history of the 
Universe. In addition, in certain wavebands AGN activity 
might be obscured b y gas and dust surrounding black hole 
(Hopk ins et al]|2006T ). More detailed consideration of these 
issues will form the basis for future work. 



ACKNOWLEDGEMENTS 

We are grateful to Zuhui Fan, Gabriella De Lucia, Roderik 
Overzier and Qi Guo for their detailed comments and sug- 
gestions on our paper. Lan Wang would like to acknowledge 
the supports from NSFC under grants 10373001, 10533010, 
and 10773001, and 973 Program (No. 2007CB815401). The 
simulation used in this paper was carried out as part of 
the programme of the Virgo Consortium on the Regatta su- 
percomputer of the Computing Centre of the Max-Planck- 
Society in Garching. 

This paper has been typeset from a TfrX/ LTpX file prepared 
by the author. 



REFERENCES 

Barger A. J., Cowie L. L., Mushotzky R. F., Yang Y., Wang 
W.-H., Steffen A. T., Capak P., 2005, AJ, 129, 578 

Colless M., Dalton G., Maddox S., Sutherland W., Norberg 
P., Cole S., Bland-Hawthorn J., Bridges T., et al., 2001, 
MNRAS, 328, 1039 

Croton D. J., Springel V., White S. D. M., De Lucia C, 
Frenk C. S., Gao L., Jenkins A., Kauffmann G., Navarro 
J. F., Yoshida N., 2006, MNRAS, 365, 11 

De Lucia G., Blaizot J., 2007, MNRAS, 375, 2 

De Lucia G., Springel V., White S. D. M., Croton D., Kauff- 
mann G., 2006, MNRAS, 366, 499 

Duschl W. J., Strittmatter P. A., 2002, in Green R. F., 
Khachikian E. Y., Sanders D. B., eds, IAU Colloq. 184: 
AGN Surveys Vol. 284 of Astronomical Society of the 
Pacific Conference Series, The Formation and Feeding of 
Massive Black Holes in the Early Universe, pp 343 — V 

Fan X., Strauss M. A., Schneider D. P., Gunn J. E., Lupton 
R. H., Becker R. H., Davis M., Newman J. A., et al., 2001, 
AJ, 121, 54 

Ferrarese L., Merritt D., 2000, ApJ, 539, L9 

Gebhardt K., Bender R., Bower G., Dressier A., Faber 
S. M., Filippenko A. V., Green R., Grillmair C, Ho L. C, 
Kormendy J., Lauer T. R., Magorrian J., Pinkney J., 
Richstone D., Tremaine S., 2000, ApJ, 539, L13 

Haehnelt M. G., Rees M. J., 1993, MNRAS, 263, 168 



Hasan P., 2007, ArXiv e-prints, 707 

Hasinger G., Miyaji T., Schmidt M., 2005, A&A, 441, 417 
Heckman T. M., Kauffmann G., Brinchmann J., Chariot 

S., Tremonti C, White S. D. M., 2004, ApJ, 613, 109 
Ho L. C, Filippenko A. V., Sargent W. L. W., 1997, ApJS, 

112, 315 

Hopkins P. F., Hernquist L., Cox T. J., Di Matteo T., 

Robertson B., Springel V., 2006, ApJS, 163, 1 
Kauffmann G., Haehnelt M., 2000, MNRAS, 311, 576 
Kauffmann G., Heckman T. M., Tremonti C, Brinchmann 
J., Chariot S., White S. D. M., Ridgway S. E., Brinkmann 
J., Fukugita M., Hall P. B., Ivezic Z., Richards G. T., 
Schneider D. P., 2003, MNRAS, 346, 1055 
Kormendy J., Richstone D., 1995, ARA&A, 33, 581 
Li C, Kauffmann C, Wang L., White S. D. M., Heckman 

T. M., Jing Y. P., 2006, MNRAS, 373, 457 
Magorrian J., Tremaine S., Richstone D., Bender R., Bower 
G., Dressier A., Faber S. M., Gebhardt K., Green R., Grill- 
mair C, Kormendy J., Lauer T., 1998, AJ, 115, 2285 
Merloni A., 2004, MNRAS, 353, 1035 
Mihos J. C, Hernquist L., 1996, ApJ, 464, 641 
Mueller A., Hasinger G., 2007, ArXiv e-prints, 708 
Richards G. T., Strauss M. A., Fan X., Hall P. B., Jester 
S., Schneider D. P., Vanden Berk D. E., Stoughton C, et 
al., 2006, AJ, 131, 2766 
Spergel D. N., Verde L., Peiris H. V., Komatsu E., Nolta 
M. R., Bennett C. L., Halpern M., Hinshaw G., et al., 
2003, ApJS, 148, 175 
Springel V., White S. D. M., Jenkins A., Frenk C. S., 
Yoshida N., Gao L., Navarro J., Thacker R., et al., 2005, 
Nature, 435, 629 
Springel V., White S. D. M., Tormen G., Kauffmann G., 

2001, MNRAS, 328, 726 
Steffen A. T., Barger A. J., Cowie L. L., Mushotzky R. F., 

Yang Y., 2003, ApJ, 596, L23 
Volonteri M., Haardt F., Madau P., 2003, ApJ, 582, 559 
Wyithe J. S. B., Loeb A., 2003, ApJ, 595, 614 



